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ABSTRACT 

Determining the properties of starbursts requires spectral diagnostics of their ultraviolet ra- 
diation fields, to test whether very massive stars are present. We test several such diagnostics, 
using new models of line ratio behavior combining Cloudy, Starburst99 and up-to-date spectral 
atlases (Pauldrach et al. 2001; Hillicr & Miller 1998). For six galaxies we obtain new measure- 
ments of Hei 1.7 /im/BrlO, a difficult to measure but physically simple (and therefore reliable) 
diagnostic. We obtain new measurements of Hei 2.06 /im/Br7 in five galaxies. We find that 
Hei 2.06 /im/Br7 and [Oiii]/H/3 are generally unreliable diagnostics in starbursts. The het- 
eronuclear and homonuclear mid-infrared line ratios (notably [Neiii] 15.6 fj,m / [Neii] 12.8 /xm) 
consistently agree with each other and with Hei 1.7 /im/BrlO; this argues that the mid-infrared 
line ratios are reliable diagnostics of spectral hardness. In a sample of 27 starbursts, [Ne III] /[Nell] 
is significantly lower than model predictions for a Salpeter IMF extending to 100 Mq. Plausi- 
ble model alterations strengthen this conclusion. By contrast, the low-mass and low-met allicity 
galaxies II Zw 40 and NGC 5253 show relatively high neon line ratios, compatible with a Salpeter 
slope extending to at least ~ 40-60 M0. One solution for the low neon line ratios in the high- 
metallicity starbursts would be that they are deficient in > 40 M0 stars compared to a Salpeter 
IMF. An alternative explanation, which we prefer, is that massive stars in high-metallicity star- 
bursts spend much of their lives embedded within ultra-compact Hii regions that prevent the 
near- and mid-infrared nebular lines from forming and escaping. This hypothesis has important 
consequences for starburst modelling and interpretation. 

Subject headings: infrared: galaxies — galaxies: starburst — galaxies: stellar content — techniques: 
spectroscopic 

1. INTRODUCTION 

In the very local (D < 5/ij^qq) Universe, the circumnuclear regions of just four galaxies (M82, NGC 253, 
NGC 4945, and M83) are responsible for ^ 25% of the current massive star formation (Heckman 1997). 
In these "circumnuclear starburst galaxies" , the star formation is confined to the inner 0.2 to 2 kpc, in a 
dense, gas-rich disk where star formation rates can reach 1000 Mq yr~^ (Kennicutt 1998). If the starburst 
initial mass function (IMF) includes significant numbers of low-mass stars, then each starburst is currently 
building up the stellar component of its host galaxy as well. A starburst enriches and heats its interstellar 
medium, as well as the local intergalactic medium. Starbursts can also drive large-scale winds that eject 
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interstellar gas, presumably casting metals into the voids and heating the gas between galaxies. Starburst 
galaxies thus play a number of important roles in galaxy evolution. 

If starbursts could be dated, then a sequence could be pieced together, charting starburst evolution 
from triggering to post-starburst quiescence. Starburst ages are most directly determined by understanding 
the population of rapidly evolving massive stars. The feedback effect of a starburst on its gas supply is 
transmitted through massive stellar winds and supcrnovac driven supcrwinds. Thus, understanding the 
evolution of starbursts and their effects on the interstellar and intergalactic media both critically depend on 
understanding the populations of massive stars. 

Unfortunately, since starburst galaxies are too far away to count individual stars, the high-mass IMF 
must be determined indirectly, in ways that are model-dependent and crude. Leitherer (1998) reviews these 
techniques and divides them into three categories: techniques to determine a lower mass cutoff by measuring 
the mass to- light ratio; to find the slope of the IMF above ~ 10 Mq; and to determine an upper mass cutoff 

from the hardness of the ionizing radiation field. 

The ionizing spectrum is set by the starburst's age, IMF, and star formation history. Consequently, the 
shape of the ionizing field spectrum is an important boundary condition on starburst models (Rieke 2000). 
Ionizing continua are often parameterized by an effective temperature (Te//), as if one stellar spectral 
type were responsible for the flux. The UV spectrum cannot be measured directly because little ionizing 
continuum radiation escapes from a starburst (Leitherer et al. 1995). Instead, the presence or absence of 
massive stars must be inferred using spectral diagnostics. Extinction in circumnuclear starbursts means that 
infrared diagnostics are preferred. 

Many line ratios have been used to estimate starburst Te//: forbidden line ratios, mixed forbidden- 
recombination line ratios like [OiiiJ/H/?, optical lines of Hei, the near-infrared line Hei 2.06 /im, and mid- 
infrared fine structure lines for example. Unfortunately, these diagnostics disagree by 2, 000 to 5, 000 K (Vanzi 
& Rieke 1997; Thornley et al. 2000; Rieke 2000), and suffer variously from intrinsic faintness, susceptibility 
to shocks and reddening, dependence on nebular conditions, and uncertain atomic constants. 

In this paper, we use the ratio of Hei 1.7 /im to Brackett 10 (BrlO) to diagnose the hardness of 
starburst ionizing fields. The faintness of the Hei 1.7 /im line restricts its measurement to nearby galaxies 
with strong emission lines. In these galaxies, the Hei/BrlO ratio should allow estimates of Te// that are 
largely independent of reddening or nebular conditions. We then use Hei/BrlO to assess the accuracy of 
diagnostics that can reach distant galaxies. Using diagnostics we find reliable, we confirm that few massive 
starburst galaxies have high-excitation spectra. While this may occur because the IMF is biased against 
high-mass stars, we propose that high-excitation spectra are scarce because the massive stars spend most 
of their main sequence lifetimes embedded in ultracompact Hii regions. 

2. OBSERVATIONS, DATA REDUCTION, AND CALIBRATION 

To assist in evaluating Tg// diagnostics, we have obtained new measurements of Hei 1.7 /im/BrlO. This 
diagnostic is unaffected by metallicity, shocks, or level pumping. Regrettably, the Hei 1.7 /im line is very 
weak, less than 10% the intensity of Br7. Thus, the Hei 1.7 /im/Br 10 ratio can only be measured in nearby 
starbursts with bright lines. 

For our sample, we chose six nearby starburst galaxies with large measured Br 7 fluxes and, when 
possible, supporting observations in the literature such as mid-infrared spectra. Near-infrared spectra were 
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obtained on the nights of 2001 April 6 and 7, using the FSPEC near-infrared spectrometer (Wilhams et al. 
1993) on the Steward Observatory Bok 2.3 m telescope. 

Table 2 lists target objects and integration times. All observations were taken with the 600 lines mm^^ 
grating, which produces effective resolutions of i? « 2000 at 1.7 /xm and R « 3000 at 2.1 fim. The slit was 
2.4" by 90". All exposures were guided by hand using an H-band camera that images the mirrored slit. 
Spectra of the six targets were obtained in H band (1.7 /xm). K-band (2.1/im) spectra were also obtained 
unless high quality spectra already existed in the literature. 

The angular sizes of the nuclear starburst regions are small compared to the length of the FSPEC slit. For 
each integration, nuclear spectra were obtained at four successive positions along the slit. For the calibration 
stars, six spectra were taken along the slit. (For brevity, we will call each resulting two-dimensional spectrum 
a "frame", and each group of frames in an integration a "set".) 

2.1. Data Reduction 

The infrared Hi and Hei lines are faint, making the data reduction approach critical. We therefore 
describe it in detail. The data were reduced using iraf First, dark frames were subtracted from the 
object frames. Each frame was then flat-fielded using a median-averaged lamp flat. The airglow and bias 
were removed by differencing neighboring frames. For the first and last frames of a set, the neighbor was 
subtracted. For each middle frame, the mean of the immediately-bracketing frames was subtracted. For the 
calibration stars the subtraction was simple. Over the longer integration times required for the galaxies (~ 4 
minutes per frame), the sky background is variable. Accordingly, prior to subtraction, we scaled each galaxy 
frame by a constant, generally within a percent of unity, to optimize the sky cancellation. This technique 
of differencing neighboring two-dimensional spectra usually removes the sky emission lines accurately. For 
a few galaxies, the resultant sky subtraction was not adequate. In these exposures, the sky lines in a set of 
frames were offset along the dispersion axis by 0.005 to 0.01 pixels, suggesting a slight, nionotonic shift in 
the grating tilt. To improve the cancellation, we used onedspec. identify and onedspec. reidentif y to 
fit, for each frame, a linear shift in the position of the airglow lines with respect to a reference frame. We 
then used images . imgeom. imshift to shift the frames, using linear interpolation, to zero the offsets. This 
was done for H-band frames of NGC 4861, He 2-10, NGC 4102, and NGC 3504 as needed. 

The next step in the reduction was to combine the frames within a set. Offsets were; detcirmincd 
by summing each frame down the dispersion axis and measuring the location of the continuum peak in 
the resulting one-dimensional image. Frames were magnified by a factor of six to permit fractional pixel 
shifts, which minimizes smearing of the data and maximizes preservation of flux during the next step of 
straightening. Magnifled images were then remapped to make the spatial axis perpendicular to the dispersion 
axis. This remapping is accomplished by fitting the tilt of the dispersion axis as a cubic polynomial; this 
function is well-defined and does not vary with time. Straightening the spectra in this manner simplifies 
subsequent extraction, stacking, and wavelength calibration (Engelgracht 1996). 

Frames were then median combined with imcombine using no rejection, and using scale and weight 
factors determined from each frame's median continuum strength (generally within 20% of unity). Frames 
were then de-magnified. Bad pixels were replaced with the average of their immediate two neighbors along 
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the dispersion axis. 

A similar procedure created two-dimensional sky frames suitable for wavelength calibrations. The 
frames were straightened and median combined as before, but without weights, scales, or offsets. The 
resulting images have no continua, only sky lines. A high signal-to-noise 1-D sky spectrum was then 
created by taking, at each wavelength, the median value over all spatial positions. 

We then extracted the spectra of the galaxies and calibration stars. Using apall in iraf , we traced 
each continuum interactively using low-order Legendre polynomials. The aperture width was chosen to be 
where the signal at the brightest part of the continuum dropped to 30% of peak. No extra background 
was subtracted at this point; trials with additional subtraction (to remove residual sky lines) added more 
noise than they eliminated. Spectra were extracted from the sky frames by using the apertures fitted for the 
corresponding objects. 

Next, the spectra were corrected for telluric absorption. To do this, each galaxy spectrum was first 
divided, in pixel space, by the spectrum of a calibration star. When a target observation was bracketed by 
calibration star observations, the average of the two stellar spectra was divided into the target spectrum. 
Otherwise, the calibration spectrum taken closest in time to the target was used. Calibration stars were 
dwarfs of spectral types F6 to GO, most within 5° of the target object. Such stars have relatively feature- 
less intrinsic spectra at ~ 2 yum, so their observed spectra reflect the variable absorption of the Earth's 
atmosphere. 

We wavelength-calibrated the sky spectra of the galaxies, using onedspec . identify and tabulated 
vacuum wavelengths of the OH lines calculated by C. Kulesa (1996, private communication). Scatter in the 
wavelength calibration of the galaxies was less than ±0.2 A, and usually below ±0.1 A. 

Wavelength solutions of the sky spectra were transferred to the [target/calibrator] spectra. Due to the 

longer integration times, the sky spectra derived from galaxy frames give wavelength solutions more precise 
than those derived from sky spectra of the star frames. The wavelength solution is quite stable with time. 

The calibration stars are not completely featureless; their metal absorption lines produce spurious 
emission lines in the [target /calibrator] spectra. To correct for this effect, we multiplied the [target /calibrator] 
spectra by the solar spectrum, which has been corrected for atmospheric absorption (Livingston and Wallace 
1991). The solar spectrum was first normalized and broadened to 10 A in H and 9 A in K to match our 
observations. Since the calibration stars are similar in spectral type to the Sun, the solar multiplication 
cures the final spectrum of the metal lines and the global Rayleigh-Jeans slope that the calibration star 
introduced. The resulting spectrum is [target/calibrator] * [©] . This procedure is explained more fully by 
Maiolino, Rieke, & Rieke (1996). 

At this point, NGC 4861 and NGC 4214 still showed residual sky lines. The spectra were improved by 
subtracting a scaled, extracted sky spectrum. For both galaxies, Br 11 and Hel 1.7 /zm are uncontaminated 
by OH lines, whereas Br 10 is somewhat contaminated in NGC 4861 and seriously contaminated in NGC 4214 
(the lowest redshift galaxy in our sample). Accordingly, we will consider only Br 11 in lieu of both Br 10 and 
Br 11 for these two galaxies. 
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2.2. Combining Spectra 

For the H-band spectra, we observed targets for 0.5 hr between calibrators. Because our total integration 
times on each galaxy were substantially longer than this, the spectra must be combined. To do this, we used 
onedspec.scombine, median combining groups of > 4 images, and average combining otherwise. 

The Hci 2.06 /im and Br7 lines were covered by different grating settings. We spliced together the two 
grating settings for each target by scaling the 2.085 /xm spectrum by a constant until it matched the flux 
level of the 2.15 /xm spectrum in the region of overlap. Scaling constants were between 1.08 and 1.3. We 
then joined the spectra using scombine. 

2.3. Subtracting the Stellar Continuum 

The nebular lines we seek sit atop a stellar continuum, whose absorption lines can mask or alter the 
emission line ratios. Before faint emission line fluxes can be measured, the stellar continuum must be 
removed. 

The lines of interest at K (Br7 and Hei 2.06 /xm) are strong enough that subtraction of the stellar 

continuum is not necessary. For H band, we used a stellar continuum template made by combining the 
15 stars from a stellar atlas (observed with the same spectrometer: V. D. Ivanov, in preparation) which 
minimized the residuals when subtracted from NGC 253 (Engelbracht et al. 1998). Thus, the template was 
chosen to be a good fit to a ~solar metallicity starburst galaxy, and was not made by modelling the stellar 
populations of each galaxy in our sample. The stars in the template arc of stellar types KO to M3, with 
metallicities between solar and half-solar. Two of the stars are supergiants, five are bright giants, five are 
giants, and three are dwarfs. 

The resolution of the template is 130 km s~^ in H. To match the intrinsic velocity dispersion of the 
galaxy spectra, we convolved this template with a Gaussian kernel as necessary to lower the resolution. 
He 2-10 and NGC 3077 were best fit with no convolution; NGC 3504 required a template with 150 km s^^ 
resolution, and NGC 4102 required a 175 km s"^ template. For NGC 4214 and NGC 4861, continuum 
subtraction was unimportant because the continua are much weaker than the nebular emission lines. 

The stellar continuum template has an absorption feature at 1.7010 ^.m. By comparison with spectra 
of the Sun (Livingston and Wallace 1991) and Arcturus (Montgomery et al. 1969), we identify this feature 

as a blend of three nickel lines and four (weaker) iron lines. At a resolution of 130 km s"-"^, the blend has an 
equivalent width of < 1.5 A, and dips to 94% of the continuum level. Because the absorption feature is a 
metal line blend, its strength will depend on metallicity. 

Besides a simple subtraction of the stellar continuum, we also added a 20% featureless continuum to the 
template, renormalized, and subtracted the new continuum from the galaxy spectra. This procedure crudely 
approximates the effect of lower metallicity. These two realizations of the continuum subtraction provide 

some estimate of the associated uncertainty. 

Continuum-subtracted H-band spectra of the galaxies and the stellar template are plotted in figure 1. 
K-band spectra are plotted in figure 2. 
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2.4. Measuring Line Ratios 

To set the continuum level, we fit a low-order Chebyshev function across each spectrum, excluding 
emission lines from the fit. Each line was fit by a Gaussian to measure the line fluxes listed in table 3. For 
the noisy, non-Gaussian line profiles of NGC 4102, we directly summed flux rather than fit Gaussians. 

To measure the relative strength of Hei 1.7 /im, wc considered two ratios: Hei 1.7 /im/BrlO, and 
[Hei 1.7/xm/Brll] X [Brll/BrlO]coseB- We assumed the value [Brll/BrlO]coseB = 0.75, which is appropriate 
for n = 10^ cm-3 and = 5,000 K (Engelbracht et al. 1998; Hummer & Storey 1987). For each galaxy, 
we computed both these ratios for both realizations of the continuum subtraction (with or without the 20% 
featureless continuum), and used the mean of these four values as the Hei 1.7 /j,m/BrlO ratio, and the 
standard deviation as an estimate of the uncertainty associated with the continuum subtraction. We also 
computed the ratio Hei 2.06 fim/Bi-y. Our measured line ratios are listed in table 4, along with values of 
Hei 2.06 /im/Br7 from the literature, and a weighted mean for Hei 2.06 /im/Br7 that combines new and 
literature values. (The quoted uncertainty for the weighted mean is the error in the mean.) 

Because of the 1.7 /im stellar absorption feature, studies that do not subtract the continuum in galaxies 
with weak Hei 1.7 /xm will somewhat underestimate the Hei 1.7 yitm line strength and therefore underestimate 
Te//. To test the magnitude of this effect, in table 4 we list both continuum-subtracted and raw (un- 
subtracted) Hei 1.7 /xm/BrlO ratios. For galaxies NGC 3504 and NGC 4102, the lines are so weak relative 
to the stellar continuum that the Hei 1.7 /xm/BrlO cannot be measured without continuum subtraction. 
For the other galaxies in table 4, the raw and continuum-subtracted line ratios are very similar; for these 
galaxies (He 2-10, NGC 3077, NGC 4214, and NGC 4861), the continuum subtraction is not an important 
source of uncertainty. 

3. MODELLING LINE RATIO BEHAVIOR 

Because the emission lines used to diagnose effective temperature have different excitation energies, 
one cannot verify that a particular diagnostic works by simply testing whether it exhibits a one-to-one 
correlation with another diagnostic. Instead, one must test diagnostics in light of photoionization models 
that, given realistic hot stellar ionizing spectra, predict line ratios appropriate to idealized nebulae. One can 
then ask a) whether the observed ratios populate the line ratio space permitted by models; and b) whether 
many observed line ratios for a particular galaxy are consistent, that is, can all be produced by one set of 
physical parameters. Thus, translating a nebular line ratio to a statement about stellar content is necessarily 
model dependent. 

Past studies have run series of models in which a single main sequence star photoionizes a nebula, 
producing tabulated line ratios as a function of stellar Te//. A measured galactic line ratio is then translated 
into an effective temperature using this tabulation (Doyon et al. 1992; Achtermann & Lacy 1995; Vanzi & 
Rieke 1997; Beck et al. 1997; Forstcr Schreiber et al. 2001). This method has its uses: namely, to compare 
model inputs and assumptions, and understand what line ratios different stellar classes can produce. However, 
in § 3.2 we will argue that, especially for mid-infrared line ratios and Hei 1.7 /um/Br 10, starbursts are poorly 
approximated by single main sequence stars; to translate a line ratio into a meaningful statement about a 
stellar population, one must consider the flux from all the stars as a function of time. First, though, we will 
consider the insights and limitations of simple one-star models. 

All our models (single-star and population synthesis) are radiation-bounded thin shells created by the 
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photoionization code Cloudy 94.00 (Ferland 1997). Using parameters determined by Forster Schreiber et al. 
(2001) for M82, we set the total hydrogen number density to Uh = 300 cm~^ and inner radius to i? = 25 pc. 
This choice of radius produces line ratios within 2% of the plane-parallel (i? = oo) case. This is because the 
shell is thin. Thus, the choice of radius only slightly affects the models, which are effectively plane-parallel. 
For single star models, we use a constant ionization parameter of \ogU = —2.3. We ran two sets of models, 
one with gas-phase abundances of solar, and the other with 1/5 solar abundances ("the low-Z models"); 
neither abundance set includes depletion onto dust grains. (We address the effects of dust in § 5.4.) Because 
Cloudy does not predict the intensity of Hei 1.7 fim, we scaled the intensity from Hei 4471 A, which shares 
the same upper level. For Case B and Tg = 5,000 K, the Hei 1.7 fxm line is a factor of 7.4 x 10~^ fainter 
than Hei 4471 A. 



3.1. Models Using Individual Stars 

We first consider the ratio of [Neiii] 15.6 /im to [Neii] 12.8 /xm. For reference, it requires 22 eV to make 

singly-ionized neon, and 41 cV to make doubly-ionized neon. Wc took ionizing spectra from the O star 
models of Pauldrach, Hoffmann, & Lennon (2001), as prepared by Smith et al. (2002), and also the CoStar 
model spectra of Schaerer & de Koter (1997) as hardwired in Cloudy. These two stellar libraries predict 
dramatically different line ratios. Dwarf, giant, and supergiant Pauldrach stars all produce a maximum 
[Ne III] /[Neil] ratio of 10 at Te//= 50,000 K. By contrast, the CoStar dwarf and giant atmospheres yield 
[Neiii]/[Neii]= 40 at Te//= 50, 000 K. At Te//= 35, 000 K, the predicted line ratios disagree by an order of 
magnitude. Forster Schreiber et al. (2001) used Pauldrach atmospheres and an earlier version of Cloudy to 
make their figure 8, which our Pauldrach models reproduce. 

The other mid-infrared line ratios also show this discrepancy. CoStar models predict ten times higher 
[Siv]/[Siii] and [Siv]/[Neii] ratios than Pauldrach models for most of the 25, 000 < Teff < 50, 000 K range; 
for [Arm] /[Am] and [Ariil]/[Neli], CoStar gives 2 and 3 times higher ratios, respectively. The near-infrared 
line ratios Hei 1.7 /um/BrlO and Hei 2.06 fim/Bvy are not sensitive to the choice of stellar atlas. 

It is sobering that current O star models predict such different mid-infrared line ratio strengths. On the 
bright side, this sensitivity suggests that mid-IR line ratios may provide astrophysical tests of O star spectral 

models in simple Hii regions. Givcon et al. (2002) performed such a test in Galactic Hii regions; they find 
that Pauldrach et al. (2001) atmospheres fit the observed [Arm] /[Am] versus [Neiii]/[Neii] relation, whereas 
stellar models that assume LTE do not. In another test, Pauldrach et al. (2001) argued that their models 
successfully reproduce the observed far-UV spectra of hot stars, as opposed to other models. Finally, Smith 
et al. (2002) argue that Pauldrach et al. (2001) atmospheres should be more realistic than those of Schaerer 
& de Koter (1997) because the latter neglect line broadening, and thus underestimate line blanketing. As 
a result, CoStar atmospheres have significantly higher ionizing fluxes, especially at energies exceeding the 
He+ edge. 

Given these problems, and that the CoStar atmospheres predict much higher [Neiii]/[Neii] line ratios 
than arc observed, we will use the Pauldrach et al. (2001) atmospheres in this paper. Still, that the Pauldrach 
spectra are better does not mean they are correct; the CoStar-Pauldrach discrepancy should serve as some 
warning of the current uncertainties regarding hot star spectra — a critical input to the models. 

We also consider Wolf-Rayet (WR) stars in simple nebulae. We use WN and WC model spectra compiled 
by Smith et al. (2002), which were generated using the code of Hillicr & Miller (1998). For a given Teff, these 
model WR stars yield much lower [Neiii]/[Neii] ratios compared to Pauldrach O stars: at Tg// = 50, 000 K, 
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the difTerence is a factor of 30 for WN, and a factor of 10^ for WC stars. This is because of the very strong 
Hne blanketing found in WR stars. Since Wolf-Rayet stars can reach much hotter temperatures than main 
sequence stars, a Tg// < 140, 000 K WC star can reach [Neiii]/[Neii] ~ 10 (comparable to the ratio produced 
by a Teff = 50,000 K Pauldrach O star). Similarly, a 120,000 K WN star can reach [Neiii]/[Neii]= 100. 
Thus, given these stellar atmospheres, only a WN star can give rise to a neon ratio between 10 and 100. 

As a result, there arc mid-in£rared line ratio regimes that only Wolf-Rayet stars can populate (again 
assuming solar metallicity.) For maximum effective temperatures of T^ff^^ — 50,000 K, Te//^^ = 
120, 000 K, and Teff^'^' = 150, 000 K, for Pauldrach et al. (2001) model O stars and Hillier & Miller (1998) 
WR stars, we find the following: 

• Main sequence O stars can only produce [Arm] /[Am] < 18, whereas WN and WC stars can reach 
ratios of 40. 

• MS O stars can only produce [Siv]/[Siii]= 0.5, while WC stars can reach 1.2, and WN can roach 3. 

• MS O stars can only produce [Siv]/[Neii]= 4, whereas WC can reach 12 and WN can reach 130. 

• MS O stars and WC stars can only produce [Ariii]/[Neii]= 2.5, whereas WN stars can reach 14. 

We have just seen that the conversion from mid-infrared line ratio to Te// is very different for main 
sequence stars than for Wolf-Rayet stars. Thus, even a modest portion of WR stars within a hot stellar 

population can significantly affect the line ratios. We also conclude that in a solar-metallicity starburst, if 
the line flux ratios exceed the maximum that main sequence stars can produce, then WR stars dominate the 
ionizing flux. 



3.2. Spectral Synthesis Models 

Given the influence of WR stars, we must consider the more realistic scenario of an evolving stellar 
population as the ionizing source. We used the spectral synthesis code StarburstQQ version 4.0 (Leitherer 
et al. 1999) to create instantaneous starbursts with an initial mass function of Salpeter-slope (Salpeter 1955) 
and initial stellar masses between 1 Mq and a variable upper mass cutoff, "M„p" (M„p= 100, 75, 60, 50, 40, 
and 30 Mq.) (M„p in this paper always refers to the IMF, not the present-day mass function.) As in our 
single-star models, this version of Starburst99 uses O star model spectra from Pauldrach et al. (2001) and 
Wolf-Rayet model spectra from the code of Hillier & Miller (1998), as prepared by Smith et al. (2002). 

We created two suites of models: the first set assumed solar metallicity in Starburst99 and Cloudy, 
and the "solar metallicity, high-mass loss" option, which is recommended as the default for Starburst99. 
(The alternative "standard mass loss" option gives qualitatively similar results.) The second set of models 
used a gas-phase metallicity of 1/5 solar in Cloudy, and the "high-mass loss, Z=l/5 solar" and "uvlines = 
Magellanic" settings in Starburst99. Dust was ignored (and will be addressed in § 5.4.) Other parameters 
were set as for the single-star models. Starburst99 calculated the spectral energy distribution (SED) of the 
burst every 0.1 Myr for 10 Myr after the starburst. The ionization parameter was normalized to a maximum 
value of logC/ = —2.3, and scaled by the number of hydrogen-ionizing photons present in the SED. Given 
the SEDs as input, Cloudy calculated line ratios as a function of starburst age. 

Figure 3 plots line ratios as a function of time for these simulations. Table 5 summarizes the spectral 
synthesis models with M„p= 100 M© and compares to line ratios from the single-star models. Line ratios 
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versus time for the low-metallicity models are shown in figure 4. 

Figure 3 shows that in the first 2 Myr, the mid infrared hne ratios fall from an initial plateau. By 
2 Myr, the 03 through 05 dwarfs (T,.// > 44, 500 K) stars in the models have left the main sequence; by 
2.5 Myr, no 03 or 04 star of any luminosity class remains. Wolf-Rayet stars, together with the remaining 
main sequence stars, create a second period of relatively high line ratios from 3.5 to 5 Myr.^ While the 
line ratios predicted for the Wolf Rayet phase are lower than predicted for WR only nebulae, clearly the 
Wolf Rayet stars are important: they produce a renaissance of high line ratios after the O stars have left the 
main sequence. It does not make sense to parameterize a mid -infrared line ratio as though the flux came 
from a single main sequence star; the ensemble of stars, including the Wolf-Rayets, must be considered. 

The models of Thornley et al. (2000), which otherwise used similar input spectra and nebular parameters 
to this work, did not include Wolf-Rayet stars. As a result, [Ne III] /[Nell] drops monotonically with time in 
their figures 6 and 10, while our solar-metallicity curves (figure 3) are double-peaked. 

4. DIAGNOSTICS OF STELLAR Te// IN STARBURST GALAXIES 

4.1. Approaches to Estimating T^/f 

In general, line ratios capable of indicating Tg// also depend on metallicity, electron temperature, 
density, ionization parameter, and the morphology of the ionized regions. Without constraints on these other 
parameters, T^// can be difficult to determine (e.g., Morisset (2003)). However, in the extreme conditions 
in starbursts, we expect less range in ionization parameter and morphology than in broad samples of Hii 
regions, and starburst metallicities can be constrained by other line ratios. Therefore, it is plausible that 
useful constraints on Te// can be derived for these regions. We return to this topic in § 5.1. 

At optical wavelengths, the ratio of [Oiii] 5007 to H/3 is frequently used as a diagnostic since both lines 
are easily observed (see for example Stasinska & Leitherer (1996).) They are relatively close in wavelength, 
and Balmer ratios can be used to correct for residual extinction. However, large optical depths of interstellar 
extinction can make [Oiii]/H/3 reflect the conditions of the outer skin of starbursts only. Thus, [Oiil]/H/? 
may not indicate the average conditions throughout a highly-extincted starburst. Also, a ratio composed of 
a forbidden metal line and a hydrogen recombination line is particularly sensitive to the metallicity, electron 
temperature, and density of the nebular region. Additionally, [Oiii] can be shock excited (Raymond et al. 
1988). 

Vflchez & Pagel (1988) have proposed an all-forbidden line diagnostic tj' , which uses lines of [Oil], 
[Oiii], [Sii], and [Siil] with wavelengths from 3726 A to 9532 A. This diagnostic was initially reported to 
work well for Hll regions (Kennicutt et al. 2000), but it is sensitive to morphology and shocks (Oey et al. 
2000). Moreover, the diagnostic is poorly suited to starburst galaxies because it involves red lines that are 
seldom observed, and is extremely subject to reddening. 

More robust optical line diagnostics can be made by comparing strengths of helium recombination lines 
with recombination lines of hydrogen (Kennicutt et al. 2000; Ho, Filippenko, & Sargent 1997; Doherty et al. 



^Thc Starburst99 model for t = 3.0 Myr and Mup = 100 Mq predicts line ratios that are sharply discontinuous from ratios 
at 2.8, 2.9, 3.1, and 3.2 Myr. (The [Nein]/[Nen] spike is 25 times higher than the surrounding points.) The M„p < 100 M© 
models and sub-solar metallicity models have no spike. Though we have been unable to pinpoint the cause from the Starburst99 
output, we feel the 3 Myr spike is spurious, not a physical effect, and we have removed it from the figures. 
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1995). He I 6678 A and Hei 4471 A are attractive for this purpose because their proximity to Ha and H/3, 
respectively, reduces reddening effects. However, these diagnostics still sample only the outer skin of the 
starburst, and the helium lines are weak, as discussed further in § 4.3. 

Because they suffer less extinction, infrared spectral diagnostics probe more deeply into a starburst 
than optical ones. For example, 10 magnitudes of extinction at 5500 A corresponds to only 1.1 magnitudes 
at 2.2 fim, and 0.8 magnitudes at 10 fim (Rieke & Lcbofsky 1985). The mid infrared fine structure lines 
are the most successful tools in this spectral region to estimate Teff (Roche et al. 1991; Kunzc ct al. 1996; 
Achtermann & Lacy 1995; Thornley et al. 2000; Forster Schreiber et al. 2001). These lines are less dependent 
on electron temperature than optical forbidden lines. However, their atomic constants are not well known 
(Fcuchtgrubcr ct al. 1997, 2001; van Hoof et al. 2000; Galavis ct al. 1997), and they are still sensitive to 
mctallicity and ionization parameter (see figure 10 of Thornley et al. 2000). 

Because of the lack of strong atomic lines, attempts to use near-infrared lines to measure Teff have 
focused on recombination lines of helium and hydrogen. For stellar ionizing sources, the hardness of the 
ionizing continuum determines the volume of He+ relative to H+ (Osterbrock 1989). For T^jf > 40, 000 K, 
the He+ and H"*" regions coincide; the Stromgren radii arc approximately equal. For lower Teffi the zone of 
ionized H extends beyond the central zone of singly-ionized He (see figures 2.4 and 2.5 of Osterbrock (1989), 
and figure 1 of Shields (1993)). Thus, by measuring the relative volumes of He+ and H+ within a nebula, 
one can constrain the effective temperature of the ionizing stellar source (s). 

The line ratio of Hei 2.06 /um/Br7 has been used to estimate Teff in starbursts (Doyon, Puxley, & 
Joseph 1992; Doherty et al. 1995) and planetary nebulae (Lmiisden, Puxley, & Hoare 2001a). However, the 
strength of the Hei 2.06 /xm 2^P — > 2^S line is not determined simply by recombination cascade, but also 
by the population in the 2^P state (Shields 1993). This level is pumped from the ground state by A = 584 A 
photons in the resonance transition l^S — > 2^P (Shields 1993; Bashkin & Stoner 1975). Photoionization of 
hydrogen, dust absorption, or Doppler shifting can change the resonance efficiency and thus the occupation 
of the 2^P state. The state can be further populated by collisions from the triplet states, primarily from 2 S 
(Shields 1993). A small Hei 2.06 /Ltm/Br7 ratio should indicate a soft continuum where there are few 584 A 
photons and few helium recombinations. Otherwise, the ratio is likely to be a poor measure of starburst 
Te// due to the dependence on nebular dust content, electron temperature, and density, as well as on the 
ionizing continuum. Some of this complex behavior is seen in photoionization models (figure Id of Shields 
1993). 

The Hei 1.7 ^m/BrlO ratio was proposed as a Tg// diagnostic by Vanzi et al. (1996), and has been mea- 
sured in several starburst galaxies (Vanzi et al. 1996; Vanzi & Rieke 1997; Engelbracht, Ricke, & Rieke 1998; 
Forster Schreiber et al. 2001) and planetary nebulae (Lumsden, Puxley, & Hoare 2001b). The Hei 1.7 /im 
line and Br 10 are close in wavelength, and A1.7 is only one-sixth of Ay, making their ratio nearly 
reddening- independent and also allowing the photons to escape from relatively obscured regions. Unlike 
Hci 2.06 /im, the Hei 1.7 nm A^D 3^P° transition arises almost entirely from recombination cascade. 
The relevant levels arc triplet states, so they cannot be pumped from the ground state, because an electron 
spin fiip would be required (Bashkin & Stoner 1975). As a result, the line ratio is insensitive to nebular 
conditions, and is determined almost entirely by the relative sizes of the H and He ionization zones. 

Figure 8 of Forster Schreiber et al. (2001) plots the behavior of the Hei 1.7 /xm/Br 10 ratio as a function 
of Teff, as predicted by photoionization models for a starburst environment ionized by hot main sequence 
stars. The Hei 1.7 /im/BrlO ratio is small for Teff< 30, 000 K because there are many more photons capable 
of ionizing hydrogen (ionization potential of 13.6 eV) than neutral helium (ionization potential of 24.6 eV). 
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For Te//> 30,000 K, the ratio rapidly increases as the zone of singly ionized helium overlaps more of the 
hydrogen Stromgren sphere. The ratio then saturates for Teff> 40,000 K, as the He"*" and H+ regions 
coincide. For rie = 100 cm~^, the saturated ratio is 

Hei 1.7/BrlO = 3.60 C1.7 [n{He)/n{H)], (1) 

where n{He) /n{H) is the gas-phase abundance of helium (by number) relative to hydrogen, and the term C1.7 
expresses the weak dependence on electron temperature (Vanzi et al. 1996). For Tg = lO** K, C1.7 = 1.000; 
other values are listed in table 1. The helium abundance n{He)/n{H) increases from the primordial value of 
approximately 0.08 (Izotov & Trihn 1998; Bono et al. 2002) to 0.1 for the Milky Way. Thus, Hel 1.7 /xm/Br 10 
should saturate at a value of 0.27 to 0.38. 

The Te// diagnostics discussed above do not necessarily agree. For example, in the starburst galaxy 
He 2-10, [Olll]/H/3 and [Oi]/Ha indicate Te//> 39,000 K (Sugai & Taniguchi 1992), whereas mid-infrared 
line ratios indicate Te//< 37,000 K (Roche et al. 1991). The Hei 2.06 /xm/Br7 observed by Vanzi & Rieke 
(1997) would indicate Te//= 39,000 K using the conversion of Doyon et al. (1992). At poor signal-to- 
noise, Vanzi & Rieke (1997) measure Hei 1.7 /im/BrlO and find Teff= 36, 000 K. This few thousand Kelvin 
disagreement translates into a serious disagreement in stellar mass: a Tg// of 36,000 K corresponds to 
approximately an 08V spectral type, which from eclipsing binaries should have a mass of ~ 22 to 25 M© 
(Andersen 1991; Ostrov et al. 2000; Niemela & Bassino 1994; Gies et al. 2002); whereas a Te// of 40, 000 K 
corresponds to an 06.5V to 07V spectral type, which should have a mass of ~ 35 Mq (Gies et al. 2002; 
Niemela & Bassino 1994). 

We now test Teff diagnostics against each other in light of the stellar synthesis models detailed above. 
Because the line physics of Hei 1.7 /xm/Br 10 is simple and well-understood (see § 4), we assume this diagnos- 
tic is unbiased, and thus accurately reflects the ionizing continuum, within the limitations of measurement 
error. 



4.2. Testing Hei 2.06 /im/Bvy 

In this section we consider the galaxies for which we obtained Hei 1.7 /im/BrlO measurements, as well 
as three galaxies with Hei 1.7 /itm/BrlO measurements available in the literature: NGC 253 (Engelbracht 
et al. 1998), for which the stellar continuum was subtracted as in this work; M82 (Forster Schreiber et al. 
2001), for which representative stellar spectra were subtracted; and NGC 5253 (Vanzi & Rieke 1997), for 
which the stellar continuum is weak enough to ignore. These three galaxies, together with the six galaxies 
for which wc observed Hei 1.7 /im/BrlO, we term our expanded sample. We also take measurements of 
Hei 2.06 /Ltm/Br7 from the literature for the galaxies in the expanded sample. 

Figure 5 plots Hei 2.06 fim/Bi-f versus Hei 1.7 ^m/BrlO. NGC 3077, NGC 4861, NGC 4214, and He 2- 
10 all have Hei 1.7 /xm/BrlO ratios consistent with the saturated value of ~ 0.3, within the measurement 
errors and the expected variation of helium abundance. Thus, these starburst regions appear to contain 
massive stars (Te//> 39, 000 K if main sequence stars.) By contrast, NGC 253, NGC 4102, and the nucleus 
of M82 have Hei 1.7 /im/BrlO< 0.15, and thus are inferred to have softer ionizing continua (Te//< 37, 000 K 
if main sequence stars.) NGC 3504 and the two off-nuclear regions of M82 have line ratios intermediate to 
these extremes. 

Figure 5 illustrates that Hei 2.06 fj.m/Bi'y does not trace Hei 1.7 /xm/BrlO as the models predict. The 
nucleus of M82 demonstrates that Hei 2.06 /im may be strong while Hei 1.7 /um is weak, contrary to the 
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expected behavior (but expected if Hel 2.06 /im is pumped.) However, for most galaxies, Hei 2.06 /xm is too 
weak for the measured Hei 1.7 /xm. This is the first direct demonstration that Hei 2.06 /im/Brj is a poor 
diagnostic of Tg/ / in starburst galaxies. Radiative transfer considerations have predicted that the behavior 
of Hei 2.06 /xm should not be a simple function of Tg// (Shields 1993). Lumsden et al. (2001a) confirm 
this complex behavior for planetary nebulae, though they attempt to constrain the dependence on Tg and 
density by also considering optical Hei lines (Doherty et al. 1995). However, the data do not contradict the 
expectation that a very low Hei 2.06 iim/Bi^ ratio (below ~ 0.2) indicates that the continuum is fairly soft, 
because there would be few ionizing photons and also few resonantly scattered photons. 

We further consider the reliability of the Hei 2.06 /Um/Br7 ratio in figure 6, by comparing it to the 
mid-infrared lino ratio [Nciii] 15.6/im/[Ncii] 12.8/im. Here, too, Hci 2.06 iim/Brj is too low for a given 
[Ne III] /[Nell] (compared to model predictions) and there is no obvious correlation between the two ratios. 
An alternative interpretation of figure 6 would be that Hei 2.06 i^m/Bij is correct and [Neiii]/[Neii] is 
systematically overproduced; we feel this is unlikely because, as we will demonstrate in § 4.4, [Nelii]/[Neii] 
is underproduced in starburst galaxies with respect to the predictions of a Salpeter IMF extending to 100 Mq. 

4.3. Testing Optical Teff Indicators 

How well do optical forbidden and recombination line ratios estimate Teff in starbursts? Figures 9 and 

10 of Kennicutt et al. (2000) show that the recombination ratios Hei 5876/H/J and Hei 6678/Ha, as well as 

[Oiii]/H/3, all track Tg// well in Milky Way, LMC, and SMC Hii regions, where Teff could be determined 
by classifying all the ionizing stars. How well do these diagnostics perform in starburst galaxies? 

In figure 7, using dereddened data from Ho et al. (1997), we compare the behaviors of Hei 6678/Hq; and 
[Oiii]/H/3 in nuclear starbursts to the predictions of Starburst99/Cloudy photoionization models. Galaxies 
with [Oiii]/H/3< 0.5 generally have low Hei 6678/Ha, indicating general agreement that Te// is low in these 
galaxies. At higher line ratios, there is considerable scatter. For most of the plotted galaxies, [Oiii]/H/3 is 
systematically high for a given Hei 6678/Hq;, compared to a solar-metallicity, M„p= 100 Mq track. Lowering 
the metallicity of the model reduces but does not eliminate the disagreement between diagnostics. Only an 
extreme model (low metallicity, M„p= 30 Mq) can fit the data reasonably well. 

A possible explanation would be that [Oiii] in starbursts is often shock-excited by supernovae (Raymond 
et al. 1988), which would be a rare effect in Hii regions and thus not affect the Kennicutt et al. (2000) plots. 
In particular, [Olil]/H/? values above ~ 2.5 require sub-solar metallicity or excitation by shocks. Thus, 
figure 7 suggests that [Oiii]/H/3 is systematically high or Hei 6678/Ha is systematically low in starburst 
galaxies. 

As a further test, figure 8 plots [Oiii]/H/? versus Hei 1.7 /xm/Br 10 for our expanded sample. Overplotted 

are Starburst99/Cloudy models as in figure 7. For low values of [OlIl]/H/3, the error bars are too large to 
judge whether the two diagnostics correlate. As in figure 7, the highest [Olli]/H/3 values observed require 
sub-solar metallicity or shock excitation of [Oill]. 

Next, we examine the behavior of [Oiii]/H/3 versus [Neiii]/[Neii] in figure 9 (omitting for now NGC 5253, 

11 Zw 40, and NGC 55 because of their low metallicity.) Optical line ratios are from the literature, and ISO 
observations of [Nelil]/[Neii] are from Thornley et al. (2000). Galaxies with [Oiii]/H/3< 0.5 generally have 
line ratios consistent with the overplotted Starburst99/Cloudy models. With higher [Olll]/H/3, the scatter 
increases. Without accurately knowing the metallicity of each galaxy in figure 9, it is difficult to judge how 
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much of the scatter in [Olil]/H/3 versus [Ne III] /[Nell] is due to the sensitivity of [Olll]/H/3 to metallicity 
rather than effective temperature. 

According to the Cloudy models, metallicity alone cannot explain the line ratios of NGC 6240, IC 1623A, 
Arp 220, NGC 3690A, and NGC 7469 (and possibly NGC 972) in figure 9. Low metallicity and a upper 
mass cutoff of 30 M© could together explain all but IC 1623 A and NGC 6240. Alternatively, aperture 
mismatch, severe extinction, or shock excitation of [Olll] could be at work. NGC 972 is not strongly 
centrally concentrated in optical emission line images, so the explanation may lie in aperture mismatch: 
the optical line ratios were measured with slitwidths of a few arcseconds, while the ISO neon lines were 
measured with a 14" by 27" aperture. The remaining discrepant galaxies all have very heavily obscured star 
formation regions, and it is likely that the discrepancy arises because the optical and mid-infrared spectra 
sample distinctly different regions along the line of sight. We also note that NGC 278 has extremely low 
[Oiii]/H/3 for its measured [Neiii]/[Neii]. Higher-spatial resolution mid-infrared spectroscopy (e.g., with 
SIRTF) may resolve this discrepancy. We will delay discussion of whether [Neiii]/[Neii] is a reliable Tg// 
diagnostic until § 4.4. 

Next, we consider the optical helium and hydrogen recombination lines, which should form niorc^ accurate 
starburst Te// diagnostics than a forbidden/recombination pair like [Oill]/H/3. To reduce reddening effects, 
we select Hei lines close in wavelength to H lines. Unfortunately, the helium lines are weak: Hei 6678 
saturates at 0.014 of the strength of Ha, and Hei 4471 saturates at 0.05 of H/3. As such, in the spectral atlas 

of Ho et al. (1997), Hei 6678 was detected in only 108 of 418 galactic nuclei, and Hei 4471 in only 16 nuclei. 
The small sample indicates that Hei 4471 is only marginally detected, and we do not consider it further. 

Figures 7 and 9 have already implicated [Oili]/H/3 as an unreliable Te// indicator for [Oill]/H/3> 0.5. 
This makes it hard to gauge the reliability of Hei 6678/Ha in figure 7. Also, the sample sizes are too small 
to compare the optical recombination line ratios to [Ne III] /[Nell], Hei 1.7 /xm/BrlO, or Hei 2.06 /um/Br7 

individually. Instead, we use the latter three Te// indicators together to test how well the optical recom- 
bination line ratios correlate with Tg//. In table 6, we list galaxies with measurements of at least two 
different Tg// indicators, in order of increasing Te//, as determined from [Ne ill] /[Neil], Hei 1.7 /xm/BrlO, 
and Hei 2.06 ^m/Br7 (when < 0.2), as available. Due to measurement error and imcertainty in the relative 
calibrations of the diagnostics, the ordering is somewhat imcertain. The published plots of the Ho et al. 
(1997) spectra lack the dynamic range to assign upper limits to the undetected optical recombination lines. 
These are marked as "non det" in table 6. 

In general, table 6 shows some correlation between Hei 6678/Hq! and Te//, though with considerable 
scatter. Using Kendall's t rank correlation test on the eight galaxies with measured Hei 6678/Ha, there is 
only a 5% chance that Te// and Hei 6678/Hq! are uncorrelated. 

4.4. Testing the Mid— Infrared Fine Structure Line Ratios 

In the mid-infrared, ratios of the fine structure lines [Neiii] 15.6 /xm, [Neii] 12.8 iiui, [Ann] 8.99 /xm, 
[Am] 6.99 /um, [Siv] 10.5 /um, and [Sin] 18.7 /xm have been used to test for the presence of hot stars 
in starbursts. From space, ISO measured these lines at low spatial resolution (14" by 27" aperture for 
[Ne III] /[Neil]) (Thornley et al. 2000; Forster Schreiber et al. 2001; Kunze et al. 1996). Ground-based 
observations (Roche et al. 1991; Achtermann & Lacy 1995) provide higher spatial resolution, but only the 
[Siv], [Ann], and [Neii] transitions can be observed through the atmosphere. As a result, ground-based 
studies must use heteronuclear line ratios, which are less ideal than homonuclear ratios available from space 
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because they are much more sensitive to elemental abundances and dust depletion. 

In table 7, we collect measurements of the mid-infrared line ratios and Hci 1.7 /xm/BrlO in starburst 
galaxies, including five regions within M82. M82 provides a testing ground for the accuracy of the mid- 
infrared line ratios as Te// diagnostics; ISO measured homonuclear line ratios in the center (Forster Schreiber 
et al. 2001), and this region has been mapped at 1" resolution in [Neii], [Arm], and [Siv] (Achtermann & 
Lacy 1995), identifying the nucleus and three infrared-bright regions nearby (all regions defined in the 
footnotes to table 7.) 

Based on the heteronuclear mid-infrared line ratios and Hci 1.7 /im/BrlO, we find that region W2 and 
the nucleus of M82 both require M„p< 65 M©, region El requires M„p< 60 M©, and region Wl requires 
M„p< 50 Mq. The heteronuclear and homonuclear mid-infrared line ratios and Hel 1.7 /xm/BrlO within 
the SWS/ISO aperture require M„p< 50 M©. Models with M„p= 100, 75, 70, or 65 Mg do not produce the 
observed ratios in any of these regions. Thus, we find that the heteronuclear line ratios give consistent ages 
and upper mass cutoffs for individual regions near the center of M82, in agreement with Hei 1.7 /xm/BrlO, 
and when averaged over the SWS/ISO aperture, give answers consistent with the homonuclear line ratios. 

We further test the mid-infrared line ratios using the five other solar-metallicity galaxies listed in 
table 7. For NGC 4102 and NGC 6240, the constraints are poor, and the line ratios can be fit by M„p= 40 
to 100 Mq. In NGC 6946, Hei 2.06 /im/Br7 and [Nciii]/[Neii] disagree unless M„p< 65 M©, but as we 
cautioned in § 4.2, Hei 2.06 fiva is not a reliable diagnostic. For NGC 253, Hei 1.7 /xm/BrlO and the neon 
ratio cannot be simultaneously matched by the M„p= 100 M© model, but models with M„p< 75 Mq can fit 
the ratios. For He 2-10, the line ratios require M„p< 65 Mq, mostly because of low observed [Ariii]/[Neii]. 

Our conclusion is that in individual regions and entire starbursts, the different heteronuclear and 
homonuclear mid-infrared line ratios and Hei 1.7 /im/BrlO give consistent answers as to age and M„p. 
This agreement supports use of the mid-infrared line ratios as diagnostics of the ionizing radiation field. 
The mid-infrared lines have large equivalent widths and a range of excitation energies, making them poten- 
tially powerful diagnostics. 

5. DIAGNOSING IONIZING CONDITIONS IN STARBURSTS 

5.1. Mid— Infrared Line Ratio Dependencies 

The mid-infrared line ratios depend on several physical parameters: mctallicity, ionization parameter, 
morphology, and the strength and shape of the ionizing continuum. To be confident in applying these ratios, 
we need to disentangle these various effects. We consider each parameter in turn. 

Metallicity. As metallicity decreases, the relative high-excitation line emission increases, because 
lower-met allicity stars have harder spectra and because lower-metallicity nebulae cool less efficiently. An- 
other effect is that Wolf-Rayct stars require larger progenitor masses with decreasing mctallicity. These 
effects can be seen by comparing the low metallicity models (figure 4) with the solar-metallicity models 
(figure 3). Using StarburstQQ and Cloudy, we find that Z= 0.2 times solar models have initial mid-infrared 
line ratios that are ^ 3 times greater than solar-metallicity models; these line ratios fall more slowly with 
time than in solar-metallicity models. 

While metallicity affects the mid-infrared line ratios, metallicity can be measured and corrected for. 
Within galaxy samples that have similar measured metallicity, uncertainties in the metallicity should affect 
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the mid-infrared line ratios by factors that are much smaller than the orders-of-magnitude changes in line 
ratio values expected due to Te// (as discussed in § 3.2). 

Ionization parameter and morphology. The ionization parameter, as the ratio of the spectral 
intensity to the gas density, combines two of the fundamental parameters that determine the degree of 
ionization in a nebula. In a Galactic Hii region, the ionization parameter changes rapidly with radius 
because of the 1/R^ falloff and absorption of UV photons by the nebula (which also alters the spectral 
shape.) Morphology then determines which parts of the nebula influence others. A starburst galaxy, however, 
is much messier than an assembly of pseudo-spherical Hii regions: the ISM is generally fragmented, and 
gas parcels are ionized by many stars. For example, in M82, it appears that the interstellar medium is 
highly fractionated (e.g., Scaquist, Fraycr, & Bell 1998) and that the whole 450 pc nuclear starburst and 
individual ^ 20 pc star-forming clusters can be described by a single ionization parameter (Thornley et al. 
(2000), citing the dissertation of N. Forster Schreiber.) Thus, it seems more appropriate to model a starburst 
as though the gas and stars are thoroughly mixed (by employing a mean UV spectrum and mean ionization 
parameter), rather than as a collection of spherical clouds, each with a single ionization source. This "mixed 
gas and stars" model is achieved in practice by assuming plane-parallel geometry and a composite ionizing 
spectrum. 

Ionization parameters (U) have been estimated in several nearby starburst galaxies by measuring the 
number of Lyman continuum photons and the size of the starburst region. Thornley et al. (2000) summarize 

measurements in NGC 253, NGC 3256, and M82. which are all consistent with log (7 = —2.3. Measurements 
have also been obtained for Arp 299 (Alonso-Herrero et al. 2000); NGC 1614 (Alonso-Herrero et al. 2001); 
NGC 1808 (Krabbe et al. 1994); IC 342 (Bosker et al. 1997); NGC 6946 (Engelbracht et al. 1996); and 
NGC 3049 (Engelbracht et al. 1996). In addition. Ho et al. (1990) measure the Lyman continuum flux 
in fourteen nearby galaxies. Six of these galaxies have multiple measurements of log U, which gives some 
estimate of the (often considerable) uncertainty. 

In figure 10 we plot the ionization parameters derived from these studies. It should be noted that each 
of these U values is actually a lower limit, since we use the maximum radius of the starburst region to 
compute the ionization parameter. When the gas density was not measured, we assume rie = 300 cm~'^; the 
true ionization parameter scales as IP = /P300 ~ log(ne/300 cm""^). Figure 10 shows that the ionization 
parameter used in our simulations, logU < —2.3, is a reasonable average value given the measurements 
available for nearby starbursts. 

How sensitive are the mid-infrared line ratios to the ionization parameter? Reducing U in our models 

by a factor of 10 lowers the [Nelli] 15.6 /xm / [Neil] 12.8 /im line ratio by a factor of ~ 7. Therefore, if the 
ionization parameters of starburst galaxies vary by a factor of ~ 10 or more, this parameter could account 
for considerable spread in observed mid-infrared line ratios. However, there is no tendency for galaxies with 
small [Ne III] /[Nell] to have low ionization parameters in figure 10, indicating that U is not the dominant 
parameter determining this flux ratio. Comparing with the restricted range of U observed in starbursts, we 
conclude from the modeling in § 3.2 that Teff dominates variations in this line ratio in such regions. 

Starburst ISM morphologies are far too complex to reproduce in simulations; fortunately, parameter- 
ization of a starburst by a single, global ionization parameter and a mean ionizing spectrum is physically 
motivated, agrees with observations, and simplifies the problem sufficiently to allow modeling. 

Another test of the diagnostic usefulness of the mid-infrared line ratios is provided by studies of 
Galactic Hll regions. Martin-Hernandez et al. (2002) found that, in compact Hii regions, the line ra- 
tios [Neiii] 15.6 urn/ [Neii] 12.8 fim, [Arm] 8.99 fxm/ [Am] 6.99 fim, and [Siv] 10.5 /xm/ [Siii] 18.7 /xm 
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correlate very well with each other, suggesting their reliability. Morisset (2003) has also demonstrated the 
use of these lines to estimate Tg// and U in Galactic Hii regions, though as demonstrated by Morisset et al. 
(2003), outside constraints on ionization parameter and metallicity are usually necessary. 

5.2. The Spectrum of the Ionizing Radiation 

We now focus on using the fine structure line ratios to estimate the spectrum of the ionizing radia- 
tion in starbursts. As figure 3 illustrates, once an instantaneous burst is older than 6 Myr, [Neiii]/[Neii], 
[Siv]/[Siii], and [Siv]/[Neii] are so low (< 0.001) that the higher-ionization hue should not be detected. 
[Arm] /[Am] and [Ariii]/[Neii] fall off more slowly, but still require a dynamic range exceeding 100 to detect 
both lines in each ratio. Such very low line ratios are not seen in Thornley et al. (2000), which with 27 
[Neili] 15.6 /Ltm/[Nell] 12.8 /um measurements is the largest sample to date of mid-infrared fine structure 
lines in starburst galaxies. The lowest ratio detected by Thornley et al. (2000) is 0.05, and 5 galaxies have 
upper limits. The simplest explanation of this behavior is that massive stars continue to form at low rates 
after the peak of a starburst. 

In the Thornley ct al. (2000) sample, all but 3 galaxies have [Nciii]<[Neii]. The three outliers, with 
neon ratios from 1 to 12, are all low-mass, low-metallicity galaxies (NGC 55, NGC 5253, and II Zw 40). We 
will consider the higher-metallicity galaxies now, in the context of the solar-metallicity models, and defer 
discussion of the low-mass, low-metallicity galaxies to § 5.3. 

In figure 3, as M„p decreases, the line ratios decrease during the main sequence phase (because the 
ionizing spectrum softens), and the gap widens between the two phases of high line ratios (because fewer 
stars become Wolf-Rayets.) We now consider these models in light of the measured neon ratios of Thornley 
et al. (2000), which are overplotted in figure 11. 

In the Mup= 100 Mq model, for 46% of the first 5 Myr, the predicted [Neiii]/[Neii] exceeds the highest 
line ratio measured by Thornley et al. (2000) for a high-mass, ~solar-metallicity galaxy; thus, this model 
poorly fits the data. A much better fit is the Z — Zq, Mup= 40 Mq model. For only 6% of the first 5 Myr 
docs this model predict [Neili]/[Nell]> 1; for 65% of that time, it predicts neon line ratios within the range 
of the Thornley detections. The M„p= 40 Mq model fits markedly better than the M„p= 50 and 30 Mq 
models. Because one-quarter of the Thornley datapoints are upper limits (excluding the three low-mass, 
low metallicity galaxies), the M„p = 40 model is a better fit to the Thornley data than the above percentages 
indicate. 

One draws the same conclusion from continuous star formation models, as shown in figure 12. Such 
models with M„p= 100 and 75 Mq predict a constant neon ratio above 1, while the neon ratio for the 
M„p= 30 Mq model falls below the Thornley range. The M„p= 40 and 50 Mq models predict neon line 
ratios within the Thornley range; the M„p= 40 model comes closer to the median. 

These results are consistent with those of § 4.4, which found that the heteronuclear and homonuclear 
mid-infrared line ratios within four regions of M82 required M„p< 50 to M„p< 65 M0 (depending on the 
region), that He 2-10 required M„p< 65 Mq, and that NGC 253 required M„p< 100 Mq. Thus, [Ne ill] /[Nell] 
in the high-mass, solar-metallicity Thornley et al. (2000) galaxies, and a concordance of line ratios in M82 
and He 2 10, are all significantly lower than the predictions of a Salpeter IMF extending to 100 Mq. An 
IMF that is deficient in massive (> 40 Mq) stars could produce the observed line ratios. 
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5.3. Ionizing Conditions in Low Metallicity Starbursts 

We now discuss mid-infrared line ratios in low metallicity starbursts. Lowering the metallicity from solar 
elevates the mid-infrared line ratios, and fills in the valley between the MS and WR phases. Hei 1.7 /im/BrlO 
is completely saturated until the WR stars die, by contrast to its double-peaked behavior for solar metallicity. 
Overplotted in figure 4 are the line ratios for II Zw 40, NGC 5253, and NGC 55, in order of decreasing 
[Ne III] /[Nell] ratio from Thornley et al. (2000). All three of these galaxies have low metallicity: II Zw 40 has 
measured [0/H] = 0.20 ±0.01 (Diaz et al. 1990) and [0/H]= 0.19 ±0.04 (Garnett 1989), [S/H]= 0.12 ±0.03 
(Garnett 1989), and [Ne/H]= 0.3 (Martin-Hernandez et al. 2002), all linear and relative to solar abundance. 
NGC 5253 has measurements of [0/H]= 0.28 (Storchi-Bergmann, Kinney & Challis 1995) and [Ne/H]= 0.58 
(Martin-Hernandez et al. 2002). NGC 55 has measured [0/H]= 0.25 to 0.37 (Webster & Smith 1983). 

We consider the line ratios of these galaxies in light of the low-met allicity models. For NGC 55, the 
only mid-infrared line ratio available in the literature is [Nelll]/[Neii]; the observed value can easily be 
produced by any M„p from 30 to 100 M©. For II Zw 40, the observed [Neiii]/[Neii]= 12 cannot be achieved 
by Starburst99/Cloudy models with solar metallicity. With the low metallicity models, we find that the 
observed [Neiii]/[Neii], [Siv]/[Neii], and [Ariii]/[Neii] line ratios cannot be produced at any age unless 
Mup is greater than 40 Mq. The Hei 1.7 /im/BrlO ratio agrees that the ionizing field is rather hard, 
but is insensitive to M„p. For NGC 5253, unless the burst is < 0.5 Myr old, the measured [Siv]/[Neli] 
requires M„p> 40 Mq. This constraint is strengthened if we consider the [Nelll]/[Neii], [Ariii]/[Noii], and 
Hei 1.7 /im/BrlO ratios, which all predict ages within 3 to 5 Myr, for a broad range of M„p (40 to 100.) The 
Hei 1.7 /im/BrlO constraint is particularly insensitive to M„p. If one assumes this age range, the [Slv]/[Nell] 
ratio requires M„p> 60 M0. 

Thus, while high-mass, solar-metallicity starburst galaxies are seen to have lower [Ne III] /[Nell] ratios 
than a Salpeter IMF with M„p= 100 Mq predicts, the low-metallicity galaxies II Zw 40 and NGC 5253 have 
the high neon ratios expected if they contain very massive stars. 

5.4. Caveats and Assumptions 

How robust is the conclusion that the nebular line ratios indicate that most high mass, solar-metallicity 
starbursts have soft ionizing continua? First, we have assumed that the Thornley galaxies are generally 
of solar metallicity. If they were more metal-poor, this would raise the predicted line ratio curves, and 
thus increase the discrepancy between the predicted and observed ratios. The opposite effect (weakening 
our constraint) occurs if the Thornley galaxies have super-solar metallicity. Thornley et al. (2000) use 
the strong-line method to derive metallicities of 1.9 ± 1 Z© for 13 of their galaxies (excluding NGC 5253 
and II Zw 40.) This result is consistent with the metallicities from optical line ratios, but we prefer the 
Thornley mid-infrared estimate because it should be reddening-independent. Starburst99 is not optimized 
for such metallicities, but we use twice-solar models nonetheless to crudely estimate whether super-solar 
metallicities could void our result. For M„p= 100, 75, and 60 Mq, doubling the metallicity from solar lowers 
the [Ne III] /[Neil] line ratios and increases the duration of the WR phase by ^ 0.5 Myr, which brings the 
models closer to agreement with observations, but deepens the trough between the main sequence and WR 
phases to 100 x below the lowest Thornley detection. These models predict neon line ratios within the 
observed Thornley range for 40% of the first 6 Myr — little better than the solar-metallicity M„p= 100 
model. To summarize, while imccrtaintics remain because; metal rich stellar evolution is not well understood, 
current models indicate that the low line ratios observed in starburst galaxies are unlikely to be explained 
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away by metallicity effects. 

Another way to negate the conclusion would be for the high-mass, solar-metallicity starburst galaxies 
to have much lower ionization parameters than we assumed. For the observed neon line ratios to arise 
in starbursts with M„p= 100 Mq, the starbursts must have U about 10 times weaker than our assumed 
log Umax = —2.3. None of the 18 galaxies in figure 10 has a measured ionization parameter this low. 

In fact, because in our models the ionization parameter starts at \ogU = —2.3 and falls with the 
ionizing flux, the ionization parameter in our models is already fairly low. (For example, 5 Myr after a 
solar-metallicity, M„p= 100 Mq burst, the ionization parameter has fallen to logC/ = —3.15.) Thus, our 
ionization parameter assumptions are conservative in that they tend to predict low line ratios for a given 
Mup-, as a result, when comparing to observed line ratios, the models will be slightly biased toward finding 
high M„p. Thus, the choice of ionization parameter is not the reason we find generally low M„p in starburst 
galaxies; the models are actually biased against finding this result. 

For simplicity, we have modelled star formation as an instantaneous burst. Starburst galaxies arc of 
course more complicated. An instantaneous burst is the most conservative assumption of star formation 
history for the purpose of constraining M„p. As illustrated in figure 12, extended star formation or a 
series of bursts would elevate predicted line ratios above the instantaneous-burst case for most of the burst 
duration. As such, extended star formation would increase the discrepancy between the low ratios observed 
in starbursts and the high ratios predicted by high M„p models. 

Dust grains harden the ultraviolet ionizing continuum, as pointed out by Aannestad (1989). Thus, if 
dust competes for the ionizing photons, this elevates the line ratios, and our conclusions are strengthened. 
Figure 11 shows this effect in Starburst99/Mappings models with and without dust. These models were 
created using the Starburst99/Mappings III web interface, beta test version 3q (Kewley et al. 2003). That 
figure also shows that the two different photoionization codes Mappings and Cloudy, given the same input 
spectra and nebular conditions, predict very similar neon line ratios. This helps address the concern that 
our results depend on the reliability of photoionization codes and their input atomic constants. 

The other major assumptions in our work are the choice of stellar evolution tracks and hot stellar 
spectra. Had we used the (hard-spectrum) CoStar models, they would have increased the predicted line 
ratios and made the Thornley et al. (2000) galaxies seem even more deficient in high-mass stars. Thus, our 
use of the softer Pauldrach et al. (2001) atlas is conservative in terms of existing hot star models. However, 
our conclusions could be invalidated if real stars have much softer ionizing continua than Pauldrach et al. 
(2001). 

We note that NGC 3077, 4214, and 4861 now have well-measured, saturated Hei 1.7 /im/BrlO, but 
no published mid-infrared spectra. Mid-infrared spectra of these galaxies should further test the trends in 
nebular line behavior discussed in this paper (all of these galaxies would appear to fall into the low-mass, 
low-met allicity category). 

5.5. UV and Nebular Diagnostics in Conflict? 

The very massive stellar populations of a number of starburst galaxies have been constrained by ultra- 
violet spectroscopy. In cases where the burst age is more than ~ 5 Myr, the UV spectra cannot test for 
stars above 40 Mq because the most massive stars have already exploded as supernovae or evolved off the 
main sequence (e.g., Gonzalez Delgado et al. (1999)). A small number of starbursts have strong P Cygni 



- 19 - 



profiles indicative of a very young burst and tfie presence of very massive stars. Thus, there appears to be 
a tendency for ultraviolet spectra of stellar populations to indicate larger M„p than do the nebular lines 
(although the galaxy samples observed in the UV and mid-infrared hardly overlap). We now consider the 
cases of He 2-10 and NGC 3049; the UV spectra of both these starburst galaxies show P Cygni profiles, and 
nebular spectra are available (mid-infrared for the former galaxy, and optical for the latter.) 

He 2-10 is an extremely rare case of a starburst which has available ultraviolet spectra of adequate quality 
to search for P Cygni line profiles as well as high-quality mid- infrared line measurements. Although He 2-10 
is of low mass and metallicity globally, the abundances in its nuclear Hii regions are approximately solar 
(Kobulnicky, Kennicutt, & Pizagno 1999). Best fits to the UV spectrum require M„p> 60 Mq (Chandar et al. 
2003). From our modelling of the mid-infrared line ratios, wc find M„p< 65 M©. Thus, these observations 
permit a discrepancy between the diagnostics, but do not require one. 

For NGC 3049, mid infrared spectra arc not available, but optical and UV spectra are. This galaxy is 
of solar (or slightly higher) metallicity in the starburst regions (Guseva, Izotov, & Thuan 2000) although it 
is of low mass and luminosity, and hence probably of low global metallicity. Gonzalez Delgado et al. (2002) 
find that the P Cygni line profiles of Civ and Siiv in NGC 3049 require M„p> 60 Mq, and rule out ages 
younger than 2.5 Myr and older than 4 Myr. Further, they find that the UV diagnostics disagree with optical 
nebular lines as to whether massive stars are present; they fit the optical nebular lines by a M„p= 40 Mq, 
t= 2.5 Myr model — parameters which would not create the observed P Cygni profiles in the UV. Given these 
results, the authors question whether nebular line ratios can reliably indicate the presence of massive stars. 

We therefore re examine the nebular lino results for NGC 3049. In modelling these lines, Gonzalez 
Delgado et al. (2002) used an older version of Starburst99 that employed pure helium WR models (Schmutz, 
Leitherer, & Gruenwald 1992) and Kurucz O star atmospheres prepared by Lejeune, Cuisinier, & Buser 
(1997). An update of Starburst99 incorporating new stellar models (Pauldrach et al. (2001) and Hillier & 
Miller (1998), as packaged by Smith et al. (2002)) became available after submission of their paper. The 
authors note that these new stellar models would soften the ionizing spectrum and reduce the discrepancy 
with the UV results, but they did not make a detailed reconciliation. 

Using our models, which make use of these new stellar atmospheres, we re-examine the nebular lines of 
NGC 3049. In Gonzalez Delgado et al. (2002), log U is fixed with time, and varies with radius (log U = —2.58 

at Rrnarr = 100 pc) in a spherical model. This results in a generally stronger U than in our models, in which 
U falls with time. To compare with the results of Gonzalez Delgado et al. (2002), we ran new models with 
logU fixed at —2.3. This value for logJ7 is within the measurement uncertainties of uh and Q(H) of the 
value used by Gonzalez Delgado et al. (2002). This choice of slightly higher ionization parameter biases our 
test toward low values of M„p (and agreement with the results of Gonzalez Delgado et al. (2002)). 

We consider the age range 3 < t < 4 Myr, as required by the UV lines (Gonzalez Delgado et al. 
2002). Over this time period, Hei 5876/H/3 can be fit by 40 < M„p < 100 Mq, and [Qui] 5007/H/3 by 
40 < Mup < 100 Mq. [Nil] 6584/H/3 only requires M„p< 75, [Sii] 6716/H/3 and [On] 3727/H/3 can be fit 
by any M„p from 30 to 100, and [Sll] 6731/H/? and [Ol] 6300/11/? cannot be fit by any model. Thus, even 
using a high U model, we do not find that low M„p is required. The updated stellar models remove the 
inconsistency between the UV and nebular lines noted by Gonzalez Delgado et al. (2002) in NGC 3049. 
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6. DISCUSSION 

We have used a number of tests to show that the nebular Une ratio [Neiii] 15.6 /Ltm/[Neii] 12.8 yum is a 
robust measure of the hot stellar population in starbursts. The line ratio is virtually unaffected by extinction, 
and as a homonuclcar ratio involving a rare gas it is not subject to abundance variation or depletion onto 
dust. Where it can be compared to other reliable Teff indicators, the agreement is good. Since the mid- 
infrared neon lines vary over several orders of magnitude during a few million years of starburst evolution, 
measurements of moderate precision can give good Tg/ / constraints. 

The neon ratio indicates low Teff in all members of a reasonably large sample of massive, high- 
metallicity starburst galaxies (Thornley et al. 2000). We have shown that plausible modifications to the 
interpretive models (adding dust, lowering the metallicity, changing the ionization parameter) leave the 
basic constraint of low Tg/ / unchanged or strengthened. 

The conclusion from [Nciii] 15.6 /im/[Neii] 12.8 /xm contrasts with the evidence for massive, hot stars 
from P Cygni line profiles in the ultraviolet spectra of two galaxies, He 2-10 and NGC 3049. In the first 
case, we find that the infrared nebular lines are consistent with the hot stellar spectrum indicated in the 
UV. In the second case, no mid-infrared spectrum exists, and the optical spectrum of Gonzalez Delgado 
et al. (2002) does not conflict with the UV result. Thus, there is no overt conflict between the P Cygni lines 
in UV starburst spectra and the limits on Teff set by nebular lines. However, to account for the observed 
low-excitation nebular spectra of starbursts, galaxies like NGC 3049 must represent a very rare stage in 
starburst evolution. Can the UV wind observations, infrared nebular line results, and starburst models be 
reconciled, given this new constraint? 

Our calculations of the emission-line properties of starbursts are based on traditional synthesis mod- 
elling, as introduced by Rieke et al. (1980). Such modelling makes the assumption that newly-formed stars 
appear on the main sequence according to an assumed formation rate with masses given by an initial mass 
function. It has recently become popular to assume a Salpeter IMF, although Rieke et al. (1993) derived 
a very similar IMF ab initio to fit the starburst properties of M82. (Both of these IMFs diSer significantly 
from estimates of the local IMF, in that both have a substantially larger portion of massive stars). 

Assuming a Salpeter IMF extending to 100 Mq, we have shown that these models predict an early 
phase in starburst evolution, of duration 3 to 4 million years, when hot, massive stars should produce high- 
excitation emission lines. "Starbursts" are identified as-such up to ages of 15 to 20 Myr; thus, about 20% 
of active starbursts should be in the early < 4 Myr phase. However, the data of Thornley et al. (2000) show 
no starbursts in massive, high metallicity galaxies with the line ratios predicted for this early phase. One 
explanation for this discrepancy would be that the Salpeter IMF substantially overestimates the numbers of 
very massive stars. We have shown that the mid-infrared line ratios can be explained if the IMF cuts off at 
40 to 50 Mq. Parameterizing the IMF by a cutoff is an oversimplification; a substantial steepening of the 
IMF slope is probably a more appropriate description. One advantage of such an IMF is that it suppresses 
the production of oxygen, which can otherwise reach very high abundances in starbursts (Rieke et al. 1993). 

However, in addition to the indications from UV spectra that stars more massive than 40 M© can 
form in substantial numbers in starbursts, the Arches Cluster near the center of the Milky Way has a large 
population of ~ 100 Mq stars (Figer et al. 2002). (The mid-infrared line ratios in the Arches (Giveon et al. 
2002) are consistent with a burst of age 2-3 or 6 Myr in our models, assuming twice-solar metallicity.) None 
of these observations can confirm the standard assumption of a Salpeter IMF extending to 100 Mq, and 
the possibility of roUoff in the IMF toward very high masses needs to be considered in detail. However, the 
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Arches and the UV starburst results suggest it is unUkely that the lack of high-excitation emission lines can 
be explained entirely in terms of a substantial steepening in the IMF above 40 - 50 Mq. 

We have therefore searched for other causes for this behavior. Wc believe an explanation can be found 
in an incorrect assumption in the standard synthesis models: that the full luminous output of newly-formed 
stars escapes into surrounding diffuse gas. This assumption justifies modelling starbursts as traditional 
low density Hii regions. Instead, wc suggest that the majority of massive stars in starbursts spend a 
substantial part of their main sequence lifetimes embedded within dense, highly-cxtincted regions — similar 
to the ultracompact Hii regions of the Milky Way — and are thus invisible to optical, near-infrared, and 
mid-infrared nebular line studies. 

In the solar neighborhood, it appears that about 15% of the main sequence life of a massive star is spent 
within an ultracompact Hll region (Kurtz, Churchwell, & Wood 1994). Hanson, Luhman, & Rieke (1996) 
have detected in the near-infrared about half of a sample of radio selected ultracompact Hll regions. They 
conclude that the detected regions typically are obscured by Ay = 30 - 50. Since the undetected regions 
in their sample should be even more heavily obscured, we take a typical case to be Ay ~ 50. Thus, these 
objects would not contribute to the optical or near-infrared emission-line spectra of the Milky Way. The 
heavy extinction would even diminish the fluxes of the mid-infrared fine structure lines such as [Nell] 12.8 fxm 
and [Nelli] 15.6 /im by a magnitude or more. More importantly, the densities in many ultracompact Hll 
regions exceed the critical densities for these lines (e.g., 2 x 10^ cm~^ for [Neiii] 15.6 /xm). Thus, even in 
the solar neighborhood, the accuracy of traditional synthesis models would be improved by assuming that 
massive stars contribute their bolometric luminosity to the region for their entire main sequence lifetimes, 
but influence the usual T,.// indicators in emission line spectra for only 85% of their lives. 

The correction suggested above would be small for synthesis modelling of the solar neighborhood. 
However, if the ultracompact Hii region lifetimes were significantly greater, a substantial deviation from 

traditional synthesis models would be expected. For nuclear starbursts in massive galaxies, the external 
pressure is large, due to both the high density and high temperature of the interstellar medium. As a result, 
the ultracompact Hii regions of starbursts should be small and their expansion retarded compared with 
Galactic ones (De Free, Rodriguez. & Goss 1995; Garcia-Segura & Franco 1996). The gravitational field of 
the central star(s) should also play an important role, slowing the expansion further (Keto 2002). Thus, it 
is likely that the massive stars in nuclear starbursts spend a substantial fraction of their lifetimes embedded 
in high-extinction regions. It is even plausible that this phase is only terminated when these stars begin 
to lose mass in strong winds — the evolutionary phase seen in UV spectra of starbursts. This possibility is 
suggested by the failure, to date, to detect any nuclear starburst that appears younger than about 3 million 
years, based on either nebular line ratios or UV spectroscopy. 

Another indication supporting the UCHii hypothesis is that starburst models under-predict the ob- 
served bolometric luminosities of starbursts. Further evidence is that radio recombination lines and free-free 
continua in starbursts indicate substantially more extinction than indicated by the Brackctt lines. For ex- 
ample, Engelbracht ct al. (1998) deduced Ay ~ 50 to the ionized gas in NGC 253 and suggested that much 
of this gas lies in very compact Hii regions. 

In addition. Beck, Turner, & Kovo (2000) found a substantial population of sources in young starburst 
galaxies whose spectra rise from A = 6 cm to A = 2 cm, indicating self-absorbed (optically thick to electron 
scattering) emission. Parsec-scale estimated sizes and large ionizing fluxes (estimated at roughly 10^ to 10'* 
OB stars) suggest that these sources are highly obscured young superstar clusters — similar to ultracompact 
Hii regions, but containing many more stars. Ironically, these results apply to low-metallicity, low-luminosity 
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galaxies similar to those that emit high-excitation nebular lines. The other galaxies known to have non-AGN 
rising spectrum sources are NGC 5253 (Turner, Ho, & Beck 1998) and He 2-10 (Kobulnicky & Johnson 1999; 
Vacca, Johnson, & Conti 2002). Massive spiral starburst galaxies do not show this self-absorbed emission, 
though it is plausible that such sources exist but are are hidden by non-thermal emission from supernovae. 

There are a number of consequences for starburst modelling. Because the duration of the obscured phase 
may depend on the mass of the central star, it may be difficult to deduce an accurate IMF in starbursts using 
nebular diagnostics. The derivation of the IMF from fitting UV spectra would also be suspect, since there 
could be a mass/age dependence on the stars contributing to these spectra rather than their providing a 
snapshot of the integrated hot stellar population. In addition, by suppressing the signatures of the youngest 
stars, the UCHii stage will tend to make the duration of starbursts appear artificially short. 

7. CONCLUSION 

We have obtained 1.7 /xm and 2.1 /xm spectra for six nearby circumnuclear starburst galaxies to measure 
the Hel 1.7 /xm/BrlO and Hel 2.06 nm/Bvy line ratios. Simple recombination physics and independence 
from nebular conditions and extinction make Hei 1.7 /xm/BrlO an accurate diagnostic of the hardness of the 
ionizing continuum ("T(.// "). The lines are too weak to be used in distant galaxies, but in nearby galaxies 
can test more detectable but potentially problematic diagnostics. 

We present models for the behavior of the Tef f diagnostics discussed in this paper. SEDs were generated 
by the spectral synthesis code StarburstQQ, which uses the most current O star and Wolf-Rayet model spectra 
(Pauldrach et al. 2001; Hillier & Miller 1998). Although Teff remains a useful shorthand term, accurate 
models of the ionizing continuum must consider the entire population of hot main sequence and Wolf Rayet 
stars. Wolf-Rayet stars maintain high line ratio values after the O stars have left the main sequence. Our 
updated models may aid interpretation of IRS spectra from SIRTF. 

We then test whether the observed line ratios are consistent with these new models. By comparison with 
Hei 1.7 /xm/BrlO, we confirm that Hei 2.06 /im/Br7 is a problem-ridden diagnostic, as predicted (Shields 
1993). Hei 2.06 /im/Br7 also fails to correlate with the mid-infrared diagnostic [Nciii]/[Neii]. However, 
we point out that a low Hei 2.06 /zm/Br7 ratio may indicate a soft UV continuum. NGC 4102, in which 
Hel 1.7 /um/BrlO and Hei 2.06 fim/Bvj are both low, is an example. 

We test [Qui] 5007 A/H/3 as a Te// diagnostic in starburst galaxies. In comparison to Hei 6678/Ha, 
[Ne III] /[Nell], and Hei 1.7 /xm/BrlO, we find that [Olll]/H/3 is systematically elevated. While aperture 
mismatch may contribute to the poor correlation with [Neiii]/[Neii], the other two diagnostics were observed 
with apertures comparable to that for [Oiii]/H/?. We suggest that shock-excitation of [Oiii] by supernovae 
is the likely cause of very high [Olll]/H/3, and that this effect plus differing sensitivities to extinction may 
explain the considerable scatter and lack of correlation with the other Tg/ / diagnostics. 

We attempt to test optical He and H recombination line ratios. Sample sizes in the literature are too 
small to compare Hei 6678/Hq: or Hei 4471/H/3 to each other or to [Neiii]/[Ncii], Hei 1.7 ^m/BrlO, or 
Hei 2.06 fim/Bij. Instead, we estimate Teff in galaxies by the latter three diagnostics, and test for a 
correlation with Hel 6678/Ha. We find a correlation at 2a significance. We re-evaluate the optical nebular 
lines of NGC 3049, and find them consistent with the age and IMF inferred from the UV spectrum; this 
removes the discrepancy between the UV and nebular diagnostics noted previously using less sophisticated 
stellar atmospheres. 
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We show that the mid-infrared lines and Hei 1.7 /im/BrlO give consistent answers for well-studied 
starbursts, lending credibility to the mid-infrared lines' use as diagnostics of starburst ionizing fields. We 
also demonstrate theoretically that the behavior of these lines in starbursts should be dominated by Tg/ / 
for galaxies with similar metallicity. 

Having found the mid-infrared line ratios to be credible Te// diagnostics, we use them to address the 
conditions in starbursts. Thornley et al. (2000) found low values of the [Nciii] 15.6 /im / [Ncii] 12.8 fim ratio 
in their sample of 27 starburst galaxies. In the context of our models, this result would suggest that high- 
mass, solar-met allicity starbursts form fewer M> 40 Mq stars than a Salpeter IMF. Adding dust, lowering 
the metallicity, choosing a different stellar atlas, or choosing a more extended star formation history would 
strengthen this conclusion. 

However, another more likely possibility can account for this result. The relatively high density and 
temperature of the interstellar medium in nuclear starbursts should increase the duration of the ultracom- 
pact Hii region phase. As a result, most of the very massive stars may spend virtually their entire main 
sequence lifetimes embedded within dense, highly extincted regions, and thus will be nearly undetectable to 
conventional optical or near-to-mid- infrared spectroscopy. This situation will make it difficult to determine 
the high-mass IMF in starbursts. 

In contrast to the high-mass, solar-metallicity starbursts, in the low-mass, low-metallicity galaxies 
II Zw 40 and NGC 5253, high neon line ratios seem to require stars more massive than ~ 40-60 M©. 
This contrast can be understood if these galaxies form stars in regions where the interstellar medium is 
less effective at confining ultracompact Hll regions, or if the lifetimes of these regions are reduced at low 
metallicity. 
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Table 1: Dependence of C1.7 on Electron Temperature (from Vanzi et al. (1996)) 



Te (K) 


C1.7 


5000 


0.95 


10000 


1.000 


20000 


1.05 



Table 2: Log of Observations for Circnmmiclcar Starburst Galaxies 

Integration times (seconds) 



Source @1.7;um @2.08;um @2.15/im 



He 2-010 


3840 






NGC 3077 


5760 


960 


960 


NGC 3504 


5760 


480 


480 


NGC 4102 


6720 


960 


960 


NGC 4214 


1920 


480 


480 


NGC 4861 


5760 






NGC 6217 




960 


960 



Table 3: Vacuum Line Wavelengths 



Transition 


A (/im) 


Brll 


1.6811 


Hei 


1.7007 


Brio 


1.7367 


He I 


2.0587 


H2 


2.1218 


Br7 


2.1661 



Table 4. Measured Line Ratios of Starburst Nuclei. 



Object 


Hel 1.7 Atm/BrlO 


Hei 1.7 ^m/BrlO 


Hel 2.06 /im/Br7 


Hel 2.06 Mm/Br7 


Hel 2.06 //m/Br7 


source: 


continuum-subtracted 


raw, this worlc 


this work 


literature 


weighted mean 


He 2-10 


0.27 ±0.013 


0.23 ± 0.05 




0.52 ± 0.03 (1); 0.64 ± 0.09 (2) 


0.53 ± 0.03 


NGC 3077 


0.46 ± 0.085 


0.49 ±0.1 


0.49 ± 0.07 


0.59 ±0.01 (1) 


0.54 ±0.05 


NGC 3504 


0.23 ±0.065 


a 


0.49 ± 0.05 


0.27 ±0.05 (3)'' 


0.49 ± 0.05 


NGC 4102 


0.141 ±0.04" 


a 


0.08 ± 0.08 


< 0.12 (3); 0.20 ±0.03 (4) 


0.08 ± 0.08 


NGC 4214 


0.31 ±0.03'' 


0.31 ± O.OS"* 


0.55 ±0.05 


0.57 ±0.07 (1) 


0.56 ±0.04 


NGC 4861 


0.3 ± 0.02'* 


0.3 ± 0.02'*''= 




0.36 ± 0.03 (1) 


0.36 ± 0.03 


NGC 6217 






0.39 ± 0.04 




0.39 ± 0.04 



References- (1) Vanzi & Rieke (1997); (2) Doyon et al. (1992); (3) Doherty et al. (1995); (4) Engelbracht (1997). 

''Without stellar continuum subtraction, the Hel 1.7 fim line cannot be measured. 

''Value ignored in computation of weighted mean. Our spectrum has much higher signal-to-noise. 

'^For the Hel 1.7 fim feature, we directly summed flux rather than fit a Gaussian. 

''Brio was contaminated by an OH sky line in this galaxy, so the Hel 1.7 //m/BrlO ratio we quote is scaled from Hel 1.7 /xm/Brll, 

assuming [Brll/BrlO]caseS ~ 0.75. 

"Only one measurement was made, so the quoted error is an estimate based on the quality of the spectrum. 



Table 5. Model Predictions for Line Ratios. 



line ratio 


value for t < 5.5 Myr, 


peak value, 


Mb Star at 


WiN at 


at 




M„p= iUU M0 


A/T 1 nn Ai" 
M„p= iUU Mq 


rp rri nnn 17" 
ie// = 0(J, UUU K 


rn 1 fin nnn 
ie// = IUU, UUU K 


rn 1 fin nnn 
i-eff = IUU, UUU K 


[Arm] /[Am] 


> 0.6 


13 


14 


25 


35 


[Ne III] /[Neil] 


> 0.05 


7 


10 


70 


0.8 


[Siv]/[Siii] 


> 0.01 


0.32 


0.45 


2.7 


0.4 


[Arm] /[Neil] 


> 0.2 


2 


2.5 


12 


0.65 


[Siv]/[Neii] 


> 0.02 


2 


4 


90 


0.9 


Hci 2.06 ^m/Br7 


> 0.5 


0.8 


0.75 


0.25 


0.4 


Hei 1.7 ;um/BrlO 


> 0.17 


0.35 


0.35 


0.35 


0.35 



Columns 2 and 3 list results from Starburst99/Cloudy models. For comparison, column 4 gives peak line ratio values for a single-star 
model using a Te//= 50, 000 K main sequence star, and columns 5 and 6 list peak values for single-star nebulae with Wolf-Rayet stars. 



Table 6. Testing optical Tg// diagnostics. 



Galcixy 


Hel i.7 /im/orlU 


Hel 2.06 /xm/Br7 


[iNeiiiJ/[NeiiJ 


Hel 6678/llQ 


TT„ T A A'~7'\ /TT a 

Hel 4471/11/5 


NGC 3628 




< 0.11 (1) 




not det 


not det 


M82 nucleus 


< 0.16 (2) 




0.16 ± 0.04 (2) 


0.35 


not det 


NGC 6946 




0.13 ± 0.06 (1) 


0.10 (3) 


not det 


not det 


NGC 972 






0.15 (3) 


0.4 


not det 


NGC 660 




0.16 ±0.02 (1) 




0.5 


not det 


NGC 4102 


0.143 ±0.04 


0.08 ± 0.08 




0.51 


not det 


NGC 6240 




0.20 ±0.03 (1) 


< 0.39 (3) 


not obs 


not obs 


NGC 3690 






0.3 to 0.7^ (3) 


0.28 


not det 


NGC 278 






0.68 (3) 


not det 


not det 


NGC 3504 


0.23 ± 0.06 






0.36 


not det 


NGC 4214 


0.31 ±0.03 






0.83 


0.36 


NGC 3077 


0.46 ± 0.084 






0.58 


0.3 



References: (1) Engelbracht (1997); (2) Forster Schreiber et al. (2001); (3) Thornley et al. (2000). All Hel 6678/Ha and Hei 4471/H/3 flux ratios are 
from Ho et al. (1997), dereddened by us, given as fractions of the saturated values (0.014 and 0.05, respectively). Uncited near-infrared line ratios: 
this paper. "Not obs" means galaxy was not observed by Ho et al. (1997). "Not det" means galaxy was observed, but line was not detected. The 
Hel 2.06 /im/Br7 ratio is only listed if < 0.2, and therefore useful. 



Two pointings. 



Table 7. Comparison of Mid-IR and Near-IR Tg// diagnostics. 



Object 


Hei 1.7 /^m/BrlO 


Hel 2.06 ^im/Br7 


ISO Mid-IR 


Ground Mid-IR 


M82-nucleus 


< 0.13 (1) 


0.48 ±0.07 (1) 






[Ariii]/[Noii]= 0.054 ± 0.02 (2a) 
[Siv]/[Noii]= 0.014 ± 0.007 (2a) 


M82-E1 










[ArIIl]/[NoIl]= 0.06 ± 0.01 (2) 
[Siv]/[Nen]= 0.026 ± 0.008 (2) 


M82-W1/B2 


0.22 ± 0.06 (1) 


0.52 ± 0.02 (1) 






[ArIIl]/[NoIl]= 0.039 ± 0.008 (2) 
[Siv]/[Noii]= 0.015 ± 0.005 (2) 


M82-W2/B1 


0.22 ± 0.05 (1) 


0.55 ± 0.02 (1) 






[Ariii]/[Noii]= 0.07 ± 0.01 (2) 
[Siv]/[Noii]= 0.013 ± 0.004 (2) 


M82-SWS 14" X 20" 


0.2 ± 0.09 (la) 


0.51 ± 0.03 (la) 


[Noin]/[Ncii] 
[Arm] /[Am] 
[Siv]/[Sni] = 


= 0.16 ± 0.04 (1) 
= 0.26 ± 0.08 (1) 
0.11 ± 0.04 (1) 


[Ariii]/[Neii]= 0.025 ± 0.004 (2b) 
[Siv]/[Neii]= 0.007 ± 0.0015 (2b) 


NGC 253 


< 0.15 (3) 


0.36 ±0.05 (4) 


[Nem]/[Nen] 


= 0.06 (5) 


[Ariii]/[Noii]< 0.3 (6) 
[Siv]/[Noii]< 0.3 (6) 


NGC 4102 


0.14 ±0.04 


0.08 ± 0.08 






[Arni]/[Noii]< 0.5 (6) 
[Siv]/[Neii]< 0.2 (6) 


NGC 6946 




0.13 ± 0.06 (4) 


[Nc in] /[Noil] 


= 0.10 (5) 




NGC 6240 




< 0.20 ± 0.03 (4) 


[Ne III] /[Neil] 


< 0.39 (5) 




He 2-010 


0.273 ± 0.015 


0.53 ±0.03 (11,8) 






[Siv]/[Noii]= 0.03 ± 0.02 (12) 
[Ariii]/[Noii]= 0.08 ± 0.03 (12) 


II Zw40 


0.43 ± 0.05 (7) 


0.31 ± 0.02 (7) 


[Ne in] /[Neil] 


= 12.0 (5) 


[S IV] /[No II] > 7 (6) 
[Ariii]/[Noii]> 3 (6) 


NGC 5253 


0.19 ±0.03 (8) 


0.48 ±0.01 (9) 


[Ne III] /[Neil] 


= 3.5 (5) 


[S IV] /[No II] = 4 (10) 
[Ariii]/[Noii]= 0.7 (10) 



References: (1) Forstor Scliroibcr ot al. (2001); (2) Achtormann & Lacy (1995): (3) Engolbracht ot al. (1998); (4) Engolbracht (1997); (5) Thornlcy ot al. (2000); (6) estimated 
from the spectra of Rociic ot al. (1991); (7) Vanzi ot al. (1996); (8) Vanzi fc Ricko (1997); (9) Lumsdcn ct al. (1994); (10) Crowthor ot al. (1999); (11) (Doyon ot al. 1992); 
(12) Beck ot al. (1997); uncitod values arc from this paper. Positions in M82 arc as follows: nucleus is defined at 9^^51"43.6^ 69?55.'00"(1950) by (1) and (2); El is defined 
by (2) as the [Noll] knol at 9^'51"'44!, 69?55.'02" (1950); Wl is defined by (2) as the [Ncll] knot at 9'? 51" 42. 5":, +69?54.'59" (1950), which corresponds to the Br7 knot 
named B2 by (1): W2 is the [Noll] knot found al 9'.'51"'41.5!, 69?54.'57" (1950), which corresponds to the Br7 knol named Bl by (1): and SWS is the 14" X 20" aperture of 
SWS/ISO, centered on 9l'51"42.2!, 69?55.'00.7'.'. Values labelled (la) arc the average over the IG'/ X lO'/ 3D field of (1), which is similar in coverage to the SWS field. Values 
labeled (2a) are estimated from figure 7 of (2). Values labelled (2b) are the average of knots Wl and W2 in (2), which provides similar coverage to the SWS field. 
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He 2-10 NGC 3077 




1.68 1.70 1.72 1.74 1.68 1.70 1.72 1.74 

wovelength (microns) wovelength (microns) 



Stellar Continuum Template 
, 2r — ■ ■ — ■ ■ ■ ■ — ■ ■ ■ ' ■ 




= 0.8 - 
0.6 - 

0.4 [ . ^ . . . ^ . . . , . . . ^ . : 

1.68 1.70 1.72 1.74 

wavelength (microns) 

Fig. 1. — H-band spectra of the starburst galaxies in our sample. For presentation, the median continuum 
levels were set to zero, and the spectra normalized so that the maximum values were unity. The galaxy 
spectra are plotted versus observed (redshifted) wavelength. The expected positions of Br 11, Hei 1.7 /im., 
and Brio are marked. Dashed lines show the continuum fits. (Note that Br 10 in NGC 4214 is contaminated 
by a telluric OH line.) For reference, the rest-frame stellar continuum template is also plotted, without 
renormalization. 
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Fig. 2. — K-band spectra of the starburst galaxies in our sample. Spectra are normalized as in figure 1, 
and wavelengths arc observed. The expected positions of Hei 2.06 /im, H2 2.12 ^m, and Br7 are marked. 
Dashed lines show continuum fits. For NGC 4102, two different continuum fits were used. 
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Fig. 3. — Line ratios versus time for Cloudy photoionization models with Starburst99 SEDs. Models have 
solar metallicity, and assume an instantaneous burst of star formation and an IMF with Salpeter slope and 
stars with masses between 1 Mq and an upper mass cutoff (M„p). 
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Time after starburst (Myr) 





Mup=100 




Mup-75 . 




Mup-50 




Mup=40 




Mup=30 




■s, - IIZw40 




\\ NGC5253 








■-. ^ NGC55 



Time after starburst (Myr) 







Mlp=100 






Mup=75 






Mup-50 






Mup=40 






Mup=30 










1 






Time after starburst (Myr) 



Time after starburst (Myr) 



Mlp=100 - 
Mup-75 - 
Mup=50 - 
Mup=40 ■ 
Mup=30 - 



0.01 I : : 1 I 

012345678 
Time after starburst (Myr) 

Fig. 4. — Line ratios versus time for low mctallicity Cloudy/ Star burst99 models. Models assume Z= 0.2 Z©, 
an instantaneous burst of star formation, and an IMF with Salpeter slope and stars with masses between 
1 Mq and Mup. Horizontal lines indicate observed line ratios for low-met allicity galaxies II Zw 40, NGC 5253, 
and NGC 55, as given in table 7. 
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M 82 



NGois: 





,NGC|4214 



GC 48^ 1 




„ NGC ^ 102 



0.1 



0.2 0.3 
Hel 1.7 |.im/Br 10 



0.4 



0.5 



Fig. 5. — Comparison of the near-infrared line ratios. Included are the galaxies we observed as well as 
NGC 5253 (Vanzi & Rieke 1997; Lumsden et al. 1994), NGC 253 (Engelbracht et al. 1998), and three regions 
of M82 (Forster Schreiber et al. 2001). The M82 regions are: the nucleus, labeled "M82 N"; the Br7 source 
"Bl" located 10" southwest of the nucleus (the upper unlabeled point in this plot); the Br7 source "B2" 
located 5" southwest of the nucleus (the lower unlabeled point.) Coordinates for the M82 regions are given 
in the footnotes to tabic 7. Filled squares show line ratio values from Starburst99/Cloudy models with solar 
metallicity and M„p= 100 Mq, run every 0.1 Myr after an instantaneous burst. The open squares represent 
models with M„p= 30 Mq. Lowering the metallicity changes the tracks insignificantly. 



-38- 



0.001 



0.01 0.1 1 

[Nelll] 15.6 nm/ [Nell] 12.8 iim 



10 



Fig. 6. — Comparison of [Nc ill] /[Neil] with Hei 2.06 /zni/Br7. Mid-infrared ratios are from Thornley 
et al. (2000), and Hel 2.06 /^m/Br7 measurements from the literature (Doherty et al. 1995; Doyon et al. 
1992; Engelbracht 1997; Forster Schreiber et al. 2001; Genzel et al. 1995; Lester et al. 1990; Lumsden et al. 
1994; Schinncrcr et al. 1997; Vanzi et al. 1996). Starburst99/Cloiidy models are ovcrplottcd: filled squares 
show solar-met allicity models with M„j,= 100 Mq; the dashed line shows solar-metallicity models with 
M„p= 30 Mq; the dot-dashed line shows the low-metallicity, M„p= 100 M© models; and the finely dotted 
line shows the low-metallicity, M„p= 30 Mq models. 
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Fig. 7. — Testing optical indicators of Tg/ / . Circles represent the subsample of galaxy nuclei that Ho et al. 
(1997) classify as being "Hii region-like". Line fluxes are from Ho et al. (1997), dereddened by us following 

the extinction law of Cardelli, Clayton, & Mathis (1989), with Ry = 3.1. Starburst99/Cloudy models are 
represented as in figure 6. Note that the [Olll]/H/3 ratios tabulated in Ho et al. (1997) are incorrectly labeled 
as dereddened, when in fact they are observed values. 
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Fig. 8. — Testing [O ill] /H/3 against Hel l.T/xm/BrlO. The filled circles arc measurements from the literature: 
dercddened [Oiii]/H/3 measurements were taken from Osterbrock & Poggc (1987); Vaccli et al. (1997); Vacca 
& Conti (1992); we also took observed [Oiii]/H/3 ratios from Ho et al. (1997) and dereddened them following 
Cardelli et al. (1989) . For M82, we use the dereddened value of Ho et al. (1997) , averaged with observed values 
from Armus, Heckman, & Miley (1989) and Wu et al. (1998). Starburst99/Cloudy models are represented 
as in figure 6. 
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Fig. 9.— [Oiii]/H/? from the literature versus [Neiii]/[Neii] from Thornley et al. (2000) (filled circles). 
Starburst99/Cloudy models are represented as in figure 6. For clarity, we omit labels for the galaxies with 
smallest line ratios. In order of increasing [Olll]/H/3 flux (or flux upper limit), they are: IC 342, M83, 
NGC 986, NGC 7552, NGC 6946, M82, NGC 3256, NGC 4945, and NGC 253. 
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Fig. 10. Ionization parameters for focal starburst galaxies, found using Lyman continuum fluxes and 
starburst radii from the literature (see § 5.1 for references.) Data from Ho et al. (1990) are plotted as circles; 
other data are plotted as squares. The vertical dashed line shows logU = —2.3, the value we have adopted 
as typical for starbursts. The vertical dotted line is the value that would depress [Neiii]/[Neii] by a factor 
of seven relative to our models, which is what is required to make the observed line ratios consistent with a 
M„p= 100 Mq Salpeter IMF. 
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Fig. 11. — Burst models of [Neiii]/[Neii] for solar metallicity, and observed values. The curves are Cloudy 
models with varying M„p; squares show Starburst99/Mappings M„p= 100 M0 models without dust; circles 
show dusty Starburst99/Mappings models with Mup== 100 Mq. The M„p= 100 M© Cloudy and Map- 
pings models use the same ionizing spectra, but different photoionization codes and treatment of the mid- 
infrared lines; the two photoionization codes are in generally good agreement. Overplotted are the observed 
[Ne III] /[Neil] ratios of Thornley et al. (2000). The neon line ratios of the three low-metallicity galaxies 
II Zw 40, NGC 5253, and NGC 55 are not plotted; they should be considered in light of the low-metallicity 
models of figure 4, and are overplotted there. 
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Fig. 12. — Continuous star-formation models of [Ne III] /[Nell] for solar metallicity. Five differerent M„p 
values arc plotted. For reference, wc overplot (finely dotted line) the solar-metallicity, M„p= 100 Mq, 
instantaneous burst model from figure 3. 



